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We give the explicit structure of the functional governing the dynamical density and current 
fluctuations for a mesoscopic system in a nonequilibrium steady state. Its canonical form determines 
a generalised Onsager-Machlup theory. We assume that the system is described as a Markov jump 
process satisfying a local detailed balance condition such as typical for stochastic lattice gases and 
for chemical networks. We identify the entropy current and the traffic between the mesoscopic states 
as extra terms in the fluctuation functional with respect to the equilibrium dynamics. The density 
and current fluctuations are coupled in general, except close to equilibrium where their decoupling 
explains the validity of entropy production principles. 
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Fluctuation theory is at the heart of statistical me- 
chanics. About a century ago appeared the first fluc- 
tuation formulae regarding equilibrium systems, in par- 
ticular from Boltzmann's statistical interpretation of the 
thermodynamic entropy. As example, the equilibrium 
density fluctuations of a gas in large volume V at in- 
verse temperature f3 and chemical potential /i satisfy the 
asymptotic law 



-fiV [0(p,n)-0(p)] 



(1) 



where ri(^) is the grand potential and ri(/i, n) = F{n) — 
fin with F{n) the free energy, is the corresponding vari- 
ational functional; at least away from the phase coex- 
istence regime where droplet formation or nucleation 
mechanisms become responsible for a slower, surface- 
exponential decay. Yet in all cases, there appears an 
important relation between the structure of equilibrium 
fluctuations and the thermodynamics of the system, mak- 
ing the equilibrium domain exceptionally well under- 
stood. In particular, the variational principles charac- 
terising equilibrium can be understood as an immediate 
consequence of its fluctuation theory and response re- 
lations can be derived from expanding ([T]) around the 
equilibrium density tlq. 

In order to include dynamics in the fluctuation the- 
ory, Onsager and Machlup derived the generic structure 
of small time-dependent equilibrium fluctuations and ex- 
plained how their dynamics relates to the return to equi- 
librium, [ij. The ensuing linear response theory for- 
malised the general relation between equilibrium current 
fluctuations and the response in driven systems in a first- 
order perturbation theory around equilibrium. To go be- 
yond and challenged by e.g. the fast progress in nonequi- 
librium experiments on nanoscale, one soon realises a 
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lack of general principles. Moreover it would be too 
optimistic to think nonequilibrium physics based solely 
on quantities typical to equilibrium descriptions supple- 
mented with the corresponding currents. Deeply related 
to that is the lack of generally valid variational principles 
for nonequilibrium steady states, beyond the approxi- 
mate ones of minimum/maximum entropy production. 
Yet, more recently there has also been great progress. 
One well-known approach to dynamical (and especially 
current) fluctuations in open systems adds to the mod- 
els fields representing the various reservoirs that count 
the long-time statistics of associated 'charges' by Master 
equation or stochastic path methods, see e.g. [1, [li 0]- 
The hydrodynamic fluctuations for some stochastic lat- 
tice gas models have been studied in e.g. For some 
standard lattice gas models the large deviations can in 
fact be explicitly calculated, see the review Q. Up to 
now, special emphasis was put on the fluctuations of the 
current, also because of relations with a celebrated fluc- 
tuation symmetry of the entropy production, 8, 9]. 

In the present letter we come back to the basic question 
whether there is at all any systematics in the fluctuations 
beyond equilibrium or close-to-equilibrium. Can one de- 
velop a formalism that would — similarly to the equilib- 
rium scheme — establish a link between the dynamical 
fluctuations and mean (thermo-)dynamical properties of 
a system, possibly with the entropy production playing 
a role similar to the entropy at equilibrium? And could 
that also explain the appearance and limitations of the 
entropy-production variational principles on a fluctua- 
tion basis? As we have shown before, (To| . the mini- 
mum entropy production principle close to equilibrium 
follows from the fluctuation theory for the occupation 
(or residence) times, which are the relative times spent 
at different states of the system. This supports both 
the relevance of dynamical fluctuation theory for under- 
standing the status and the validity of various nonequilib- 
rium variational principles, and the importance of time- 
symmetric observables in these considerations. Indeed 
our results here strongly suggest that only by treating 
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jointly the time-symmetric and -antisymmetric sectors 
do model-specific results make place for a unique fluctu- 
ation structure. As far as we know, that is one of the 
rare occasions where the nonequilibrium world can be 
seen submitted to general laws. 

To address the above questions and in the context of a 
stochastic network we propose to study the joint dynami- 
cal fluctuations of the occupation times (time-symmetric) 
and currents (time-antisymmetric). We will show that 
these joint fluctuations have an explicit and general struc- 
ture, with a fluctuation functional derived from the so- 
called traffic measuring the mean dynamical activity in 
the system, which is hence the counterpart of the entropy 
or grand potential in the equilibrium static fluctuation 
theory, cf. ([T]). Only close to equilibrium there emerges 
a simple relation between that traffic and the entropy 
production. Together with a decoupling between small 
time-symmetric and time-antisymmetric fluctuations in 
the close-to-equilibrium domain, this lies behind the ap- 
proximate validity of the entropy production principles. 
This substantially extends the argument in [10l|. Our 
main results are relations ITTI) . PT|) . and (|77|) below. The 
application of our formalism but for driven diffusions can 
be found in QJ^. Our emphasis here on jump processes 
makes the analysis also suitable to the statistics of quan- 
tum transport as e.g. in [3, particular examples will 
follow separately. 

In contrast with full counting statistics methods, the 
reservoirs are not made explicit in our approach. Instead, 
we assume that the changes in all reservoirs or loads are 
mutually distinguishable and can be read off from the 
trajectory of the system. Another remark concerns the 
meaning of the occupation times. This is a dynamical 
observable and its fluctuations fundamentally differ from 
static fluctuations, which count the plausibility that the 
system obeys some statistics given that the system was 
in its typical stationary state far in the past. These static 
fluctuation functionals are nonequilibrium variants of the 
equilibrium free energy and have been extensively studied 
in the context of lattice gases in the hydrodynamic limit, 
[13. As a matter of principle, one expects that one could 
recover the static from the dynamical fluctuations that 
are studied here. 

The mathematics involved is the theory of large devi- 
ations and stems from the work of Donsker and Varad- 
han, [l^l . A useful and repeatedly exploited technique in 
this approach is to compute the fluctuation functionals 
on more coarse-grained levels from constrained minimi- 
sations of a fine-grained functional. That is called the 
contraction principle. 



I. GENERAL FORMALISM 

We consider a mesoscopic nonequilibrium system mod- 
elled by a stationary ergodic process Xt making jumps 
on a discrete set of states, {x, y, . . .}. As is typical for a 
thermodynamic formalism it is not essential for the math- 



ematical structure whether the process represents a single 
particle random walk or a many-body open system. We 
only ask the dynamics as given by some transition rates 
w{x, y) on ordered pairs a; — > y to be ergodic. For an easy 
interpretation we assume that the local detailed balance 
principle applies, according to which log [w{x^ y)/w{y, x)] 
is the entropy change in the environment (possibly made 
of several distinct reservoirs) per single event x y. 
Tracing the whole trajectory of the system, all currents 
as well as the total entropy exchange with the environ- 
ment can be determined. 

We start from that fine-grained level of description 
and we consider as dynamical observables the occupa- 
tion times 

Pt(x) = ^^ X{u;t^x)dt (2) 

(with X equal to one or zero, indicating whether the event 
in brackets occurs, respectively does not occur) jointly 
with the two-point correlations, for all x ^ y, 

1 

CT{x,y)St ^ — J xi^^t ^ x) x{t^t+st ^ y) dt (3) 

counting the number of jumps x y, both defined 
for each realisation of the process {uJt', < t < T). 
The occupation times pt{x) form a random distribution 
that asymptotically approaches the stationary distribu- 
tion, \miT— > 00 Pt{x) = pix), with probability one by the 
ergodic theorem. Similarly, the empirical correlations 
have the almost-sure asymptotics ImiT^oc C'Tix,y) = 
p(x)w{x,y). 

The question about dynamical fluctuations concerns 
the long-time asymptotics of possible deviations of px 
and C't from their typical values: to compute the prob- 
ability Pt(p, k) to observe for all x and y, 

Pt{x)=p{x), CT{x,y) = p{x)k{x,y) (4) 

We must add here the stationarity condition 
J2y[piy)k{y,x) — p{x)k{x,y)] = since, by conservation 
of probability, limT^oo J^yiCrix^y) - CT{y,x)] = for 
every realisation of the process. To determine Pt(p, ^) 
of ([4]), we compare the path-distribution of the original 
stationary process with rates w to a fictitious stationary 
process with rates k. The former distribution reads, 
with X{x) = w(x, y) the escape rates, 

PT(t^) = pixo) e-^^^'>^'^wixo,xi)dti e-^(^i)(*2-*i) . . . 

...u;(x„_i,a;„)di„e-^(^")(^-*") (5) 

on realisations lu — (xq, 0; xi, ti; . . . ; a;„, i„ < T) with 
jumps Xk~i — > Xk at times t^. The probability P^(uj) 
of the same realisation under the fictitious process is ob- 
tained by replacing w with k, and A with the escape rates 
J2y k{x,y). We exploit that (i) for any trajectory u sat- 
isfying the constraints ^ the density of with respect 
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to the only depends on the time-averages p{x) and 
k{x, y); (ii) those values p and k become typical under the 
fictitious process Pf. when T ^ oo. Using both proper- 
ties, the probability under study is 



PT(p,fc) = PT(P,fc) 



dP* 



conditions ^ 



dP* 



(6) 



asymptotically for T ^ oo. Explicitly, Pt(p, A:) 
exp[— TX^j, k)] with the fluctuation functional 



^{P,k) ^^p{x) k{x,y) log 



k{x,y) 
w{x,y) 



■k{x,y) + w{x,y) 



(7) 

(remember that /(p, k) — oo whenever p is not stationary 
with respect to the transition rates k). This result is our 
starting point towards a systematic generation of various 
other fluctuation laws by contraction, in both the time- 
symmetric and the time-antisymmetric domains. 



II. OCCUPATION-CURRENT FLUCTUATIONS 

The observed time-averaged currents correspond to 
the antisymmetric part of the two-point correlations, 
Crix^y) — CTiy,x). The joint fluctuation law for the 
currents and the occupation times 



-TI(p,j) 



(8) 



can be derived from ([7]) by solving the minimisation prob- 
lem 

I{P,3) = i^}i{^{p,k)\p{x)k{x,y) -p{y)k{y,x) ^ j{x,y)} 

k 

(9) 

For stationary currents, j{x, y) = 0, to which we from 
now on solely restrict, the solution I{p,j) = I{p,k*) is 
determined from 

k*{x,y)^w{x,y) e^(^'2')/2 
A{x,y)^-^{y,x) (10) 
3{x,y) ^p{x)k*{x,y) -p{y)k*{y,x) 

(Otherwise I{p,j) = oo.) As a result, 

HpJ) = " l^Ki^^y) - tp{x,y)] 

(11) 



x,y 



in which 



tpix, y) = p{x)w{x, y) + p{y)w{y, x) (12) 



and 



Kix.y) ^p(x)k*{x,y) +p(y)k*(y,x) (13) 



measure the mean dynamical activities; we call them traf- 
fic and they yield the symmetric counterpart to the ex- 
pected currents. The second term in (fTT|) is therefore an 
excess in the overall traffic needed to create the fluctua- 
tion or to make it typical. Similarly, by the local detailed 
balance principle, the first term corresponds to an excess 
in the entropy flow to the environment which amounts to 
^ = I Ea:,j/i(2;,y)log[t«(a;,?/)/w(y,x)] under the original 
process and analogously for the modified one. 

Next, being motivated by the equilibrium fluctuation 
theory, cf. ([1]), we reveal a hidden canonical structure 
that enables a particularly illuminating formulation of 
our result. Any nonequilibrium process can be related 
to a reference detailed balanced one with rates wo{x,y), 
so that w{x,y) — woixjy) e"'^^'^^^^ with some driving 
a{x,y) — —a{y,x). (For example, the rates wo{x,y) = 
\/w{x, y) w{y, x) e*^^)^'*'^) and s an arbitrary state func- 
tion, can serve as such a reference.) Having fixed wq, the 
rates w{x,y) = w^{x,y) are now parameterised by the 
driving 17(2:, y), and we introduce the potential function 

H{p,a) 2^p{x) [wa{x,y) - wo{x,y)] (14) 

equal to the excess in the overall traffic with respect to 
that reference. It is a potential for the expected transient 
currents jp,a{x,y) = p{x)wa{x,y) - p{y)w„{y,x) in the 
sense that 

5H{p, ct) = i ^ ip^a {x, y) 6(r{x, y) (15) 



(with the p kept fixed in the variation). Its Legendre 
transform is 



G{p,3) 



sup I 2 2^'^'ix,y)j{x,y) - H{p, a') 



(16) 



and we observe that the suprcmum (taken over all an- 
tisymmetric matrices) is attained at a' = a* such that 
jp^a* = J J which means that the driving a and the current 
j are canonically conjugated variables. Further, eqs. (fTUl) 
are solved with A = cr* — cr, hence the fluctuation func- 
tional I{p,j) = Ia{p,j), eq. PT|) . obtains the final form 



I4p,j) ^^[G{p,j)+H{p,a) ~ 5(a,j)] (17) 



with 

S(sf,f)^\^^a{x,y)j[x,y) (18) 

x,y 

the observed entropy current into the environment. That 
is our main result, giving the fluctuation functional en- 
tirely in terms of the potential function H(p,a) (i.e. in 
terms of the overall traffic) and derived quantities. The 
functional G{p,j) directly gives the reference equilibrium 
dynamical fluctuations as Io{p,j) = ■^G{p,j), hence (fT7)) 
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specifies the nonequilibrium correction to that equiHb- 
rium. Remark also that the antisymmetric part of the 
functional la under time reversal equals Ia-{p,^j) — 
Icr{P:j) = compare [§, [l^, [l^, which is just 

the steady state fluctuation symmetry. However, more 
important is that (fT7|) also in a generic way specifies 
the time-symmetric component. That is why (jl7p rep- 
resents a generalised Onsager-Machlup Lagrangian de- 
scribing steady fluctuations, the generalised dissipation 
functions being G and H . At the same time, one recog- 
nises the mathematical structure of equilibrium fluctua- 
tions; the grand potential and the variational func- 
tional n) of (P) get replaced here by —H{p, a)/2 and 
[G{p,j) — S{a,j)]/2, respectively. Note that while the 
potentials G and H depend on the choice of reference 
equilibrium dynamics, the resulting functional I{p, j) is 
of course independent of that. 

Fluctuation laws on a more-coarse grained level, e.g., 
the fluctuations of a single selected current, can be ob- 
tained by further contractions starting from ([7]) or ([T7]). 
Then, depending on the particular question, a modified 
canonical formalism can be established. 

There is a trivial yet important generalisation of the 
above results to systems in which a transition x — > y can 
go via multiple channels, each possibly corresponding to 
the interaction with different reservoirs. For these sys- 
tems the formulas (O, pT|) . (fTS]) . etc remain valid 
if the ordered pairs x,y in the sums get replaced with 
x,y,a, the a labelling the channels. A simple example 
of such a multi-channel model comes in the next section. 



III. EXAMPLE 

We demonstrate the above formalism on a model of 
transport over a single level (quantum dot). There are 
two configurations a; = 0, 1 corresponding to the level 
being empty respectively occupied, and it is coupled to 
the left (L) and the right (R) reservoirs. Using the nota- 
tion Vl and Vr for the potential gradients between that 
level and the reservoirs, both oriented in the L ^ di- 
rection, the local detailed balance principle restricts the 
possible transition rates corresponding to each channel 
to the following general form: 



(19) 



For simplicity, we consider here only the case = Tr = 
F. Writing the occupation times as p{0) = (1 — w)/2 
and p{l) = (1 + f)/2, the expected transient currents 
(also both oriented in the L ^ R direction) and traffic 
separately for each channel equal 



j^' _ p gjj^]^ _ -p Yv cosh ■ 



tl;-^ = Fcosh 



L,R 



=p Tv sinh 



2 

f3VLM 



(20) 
(21) 



As a reference equilibrium we take the dynamics (fTQ]) 
for Vl = Vfi — (with the symmetric part F kept un- 
changed) . The current potential ([T4|) is determined from 
the overall traffic: 



H{v,Vl,Vr) = 2T (cosh 
-I- cosh 



PVl . .I3Vl 

V smh 

2 2 

I3Vr 



V smh — 2 



(22) 



One checks that dH/dVLM. = Pjt^' which IS an in- 
stance of ([Tsl) . The Legendre transform of H at = 
— j gives the occupation-current fluctuation func- 
tional G{v,j) — /o(f,j)/2 for the reference equilibrium 
dynamics, cf. (fT6|) : 



Giv,j) 



sup [/3j{VL 

Vl,Vb 



VR)^Hiv,VL,VR)] 



= 4j log 



4F 



1 



1 - 



F2 



(23) 



F2 



This extends to the nonequilibrium dynamics by the gen- 
eralised Onsager-Machlup formula (fT7|) . E.g., in the L- 
R symmetric case Vl = Vr = V, the entropy flow is 
S = 2pVj, and the nonequilibrium fluctuation functional 
becomes 



Iv{v,j) ^ 2jlog 



1 



2F 



cosh 

2 



F2 



-pvj 



(24) 



Due to the 'particle-hole' symmetry I[—v,j) = I{v,j), 
the (marginal) current fluctuations correspond to the rate 
^v(j) = lv{0,j), which is 



-I-2F 



11' 
F2_ 

cosh 



PVj 



PV 



1 



21 

F2 



(25) 



Again by contraction, the fluctuation functional for the 
occupation times is 3y(u) — Ivi'^ij*) where j* — 
T^/l — sinh(/31//2) is the most probable value of the 
stationary current given u. As a result, 

Zviv) = 2Fcosh(^) (1 - ^JX-v-^) (26) 



IV. REGIME OF SMALL FLUCTUATIONS 

The main features of the joint occupation-current fluc- 
tuations already become manifest in the leading order 
around the nonequilibrium steady state. For our origi- 
nal dynamics with stationary distribution p, steady cur- 
rent j and steady traffic t, we write p = p{l -\- ewi). 
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j = j + eji. Standard perturbation theory applied 
to lU), up to quadratic order in e, gives as a final re- 
sult I{p,j) = I^'^\ui, ji) where 



/"'("i.ji) 



-'--nV+u, + '--{y+u,f\[x,y) (27) 

with the shorthand V^ui(a;, = [^1(2;) ± Mi(y)]/2. 
This formula clearly demonstrates how the occupation 
times and current become coupled away from equilib- 
rium. That coupling is proportional to the station- 
ary current and vanishes only in the close-to-equilibrium 
regime where j — 0{e). 

The appearance/disappearance of the occupation- 
current correlation is deeply related with the valid- 
ity/breaking of the entropy production principles. The 
expected value of the (transient) entropy production rate 
£{p) at a given distribution p is the sum of the ex- 
pected entropy flow \ Y,^,y3p{x,y)\og[w{x,y)/w{y,x)] 
and the rate of increase of the system's entropy 
~5 y 2/) l'^§?'(-^)' In the same quadratic ap- 
proximation as above but now close to equilibrium so that 
w = Wti[l + 0(e)], the entropy production rate equals 



(28) 



with 2 = 0{t). On the other hand, from (ITT)) 
the marginal distribution of the occupation times for 
j — 0(e) corresponds to the functional Z^'^^i^x) = 
\ J2x,y tiV^ui)^{x, y), and hence 



(29) 



see for more details. Hence, the stationary distribu- 
tion p is a minimiser of the entropy production rate and 
the latter governs the occupation fluctuations — this is no 
longer true beyond the close-to-equilibrium regime where 
the occupation-current correlation becomes relevant. A 
similar argument reveals a direct link between the cur- 
rent fluctuations and the maximum entropy production 
principle, ^lli] . 



sition intensities and occupation times in a steady state 
regime described by a Markov jump process, ([7]). In the 
occupation-current form (|17p , it gets a remarkable canon- 
ical structure: the (reference) equilibrium functional is 
corrected by its Legendre transform which is just a po- 
tential for the expected currents, and by the entropy 
flow. These functionals form a natural starting point 
towards the study of fluctuations for any selected col- 
lection of observables that can be expressed in terms of 
transitions/currents and occupations, via the contraction 
principle. That provides an alternative to the existing 
approaches. 

As a new and crucial quantity, unseen in close-to- 
equilibrium considerations, enters the traffic, measur- 
ing the time-symmetric dynamical activity in the sys- 
tem. This observable naturally enters beyond the lin- 
ear response theory, e.g., in determining the ratchet cur- 
rent and in the escape rate theory, [l^ . The overall 
traffic yields the current potential, and its excess together 
with an excess in the entropy flow directly determine the 
joint occupation-current fluctuations, pTjl . 

The time-symmetric and time-antisymmetric fluctua- 
tions mutually couple even for small fluctuations around 
the nonequilibrium state, (|27p. Their decoupling in lead- 
ing order around equilibrium is a fundamental reason for 
the known stationary variational principles to be approx- 
imately valid. 

For extended systems with a large number of de- 
grees of freedom, phase transitions may become visible 
through singularities of the fluctuation functionals, p^ . 
It should indeed not escape the attention that the analy- 
sis from ([TT|) to pT|) requires some strict convexity argu- 
ments and uniqueness of solutions. That is certainly one 
of the most fascinating possibilities that can be discussed 
within our general framework. 
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